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Abstract 


The paper is devoted to the derivation of a simple approximate equation for the reflection function of 
semi-infinite cloudy media for the case of small and moderate observation angles. The absorption of light in 
cloud is neglected. The accuracy of equation obtained is studied using the numerical solution of the radiative 
transfer equation. It was found that approximation derived can be used for large scattering angles (larger 
than 150°) with accuracy better than 5% if observation angles are smaller than 45°. The error is below 5% 
for observation angles smaller than 30°, incidence angles smaller than 55° and arbitrary relative azimuths of 
incident and reflected light beams. 
© 2003 Elsevier Ltd. All rights reserved. 
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1. Introduction 


The reflection function of a semi-infinite nonabsorbing medium plays an important role in the 
general radiative transfer theory [1]. In particular, its knowledge allows for the calculation of trans- 
mission, reflection, and absorption characteristics of light scattering layers of large but finite thickness 
(e.g., optical thickness t > 10 for clouds), using asymptotic radiative transfer equations [2,3]. This is 
why this function has been studied extensively by many authors [1,4—8]. In particular, a user-friendly 
and highly accurate code for the calculation of this function is available on the Internet [6]. Various 
approximate solutions for this function also has been reported [1,5,7,8]. The task of this paper is to 
derive a simple and yet accurate approximation for the reflection function in a special case of water 
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clouds in visible. In essence this is the extension of our approach reported earlier [8] on the case 
of non-nadir measurements. We propose a simple parameterization for this function and check its 
accuracy. The approach proposed can be easily extended to the case of discrete random media other 
than clouds. 


2. Theory 


The reflection function of a semi-infinite layer R(C, n, o) (€ is the cosine of the observation angle 
V, is the cosine of the incidence angle Vo, and ¢ is the relative azimuth of incident and reflected 
beams) can be written as [6] 


RENo) = RCE m) +25 Rul.) cos(mg). (1) 
m=1 
Here 
z 1 2n 
Ram =5- | RGno)de Q) 
and 
1 2n 
Rll =5- f RE oos) do. 3) 


R(,n) gives the reflection function averaged over the azimuth and R,,(¢,7) is the mth Fourier 
coefficient of the reflection function. Fourier components R,,(¢,1) for water clouds were studied in 
detail by King [9]. 

Eq. (1) is a very convenient starting point for the parameterization purposes because it decouples 
(¢,7)-dependences from the azimuthal dependence ø. It is much easier to parameterise functions 
RC n) = R(n, ©) and Rm6, n) =Rn(n, ©) than the three-parameter function R(C,y,~) [7]. 

The second term in Eq. (1) is of no importance for nadir observations, where the azimuth p plays 
no role due to the symmetry of the problem [8]. For non-nadir (but not too close to a horizon) 
observations the single light scattering gives a major contribution to the second term in Eq. (1). 
This contribution can be easily found [5]. It is given by the following equation: 


_ pO) — BO) 
ar (4) 
where 
_ 1 27 
wO)=5- | Odo (5) 
0 = arccos(—(n + /(1 — EX = n?) cos @) (6) 


and p(@) is so-called phase function. This function gives a conditional probability of light scattering 
in a fixed direction and normalized as follows: 


1 1 
3 | rO snoa = l, (7) 


where 0 is the scattering angle. 
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Replacing the second term in Eq. (1) by Eq. (4) we obtain [5]: 


m 0)— p0 
RENOJ = Rn) + ae (8) 
Note that the term 
— p@®) 
EET Á 


gives the contribution of single scattering into reflection function of a semi-infinite layer [3,10]. 

In essence, Eq. (8) states that only the azimuth dependence of singly scattered light plays a 
primary role in the overall R(ọ)-dependence. This is generally not a correct statement as shown, 
e.g., by King [9] and Melnikova et al. [7]. However, such an approach can be used with a high 
accuracy in a number of viewing geometries as specified below. 

Clearly, the azimuth dependence plays no role for the nadir observation (¢ = 1) and illumination 
(n = 1) conditions. Then we have 


RO, n) = RE 1) = RC, n) = RG 1). (10) 


This means that the second term in Eq. (8) should vanish in this case. This is indeed the case. In 
particular, we have from Eq. (6) at y = 1: 0 = n — 1, independently on the azimuth and, therefore, 
P(O) = p(@) (see Eq. (5)). _ 

The parameterization of the function R(¢,7) can be achieved as follows. First of all we note that 
it follows for isotropic scattering (p(@) = 1) [10]: 


5 eln) 
R(E, n) = RCC, n) = s, 11 
CM=REM= Fe (11) 
where the function g(¢) can be approximated by the following linear dependence [3]: 
o) = a(1 + 2¢). (12) 
where œ = 4,/3/7. Substitution of Eq. (12) into Eq. (11) gives 
z A+B(C4+n)4+C 
Reny EED tN (13) 


4(¢ +n) 


where A = Ë = 1,B =24,C =44. 
Now we assume that Eq. (13) holds also for the case p(0) ¢ 1 but coefficients A,B,C then depend 
on the phase function of a particular medium. In the first approximation only the dependence of 


A,B, and C on the asymmetry parameter 


1 


s= p(@) sin 0 cos 0 dé (14) 
0 


may be considered as it was done, e.g., by Melnikova et al. [7] for a particular case of a Henyey— 
Greenstein phase function. Note that the parameter g for water clouds, considered here, varies only 
weakly in visible (g ~ 0.85) and we can assume the same constants A, B,C for clouds with arbitrary 
microphysics. 

The approximate validity of Eq. (13) is confirmed by the numerical calculations using the 
Ambartsumian’s nonlinear integral equation for R(¢,7,@). These results (for the function D(¢,7) = 
4(¢ + n)R(¢,n)) are given in Fig. 1. Indeed we see that the dependence D(n) is very close to a 
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n 


Fig. 1. The dependence D(n) at various values of ¢ for the Cloud C1 model at 2=0.7 um. 


linear function at ¢ = const and n > 0.2. However this is not the case, at 7 < 0.2 (not shown in 
Fig. 1). We see also some deviations from the linear dependence for the glory (¢ = 1,7 > 0.9) 
scattering region. 

Note that calculations in Fig. 1 were performed for the Deirmendijan’s cloud C1 model at the 
wavelength 2=0.7 um [11]. Then it follows for the effective radius of droplets aef = (a*)/(a?)=6 um 
Here (---) means averaging on the droplet size distribution f(a) = Aa°e~°"/* (a is the radius of 
droplets, Æ is the normalization constant). 

Performing the parameterisation of data presented in Fig. 1, we obtain 


A ~ 3.944, Br-25, Cr 10.664, (15) 


which differs from the isotropic scattering case, where A ~ 1,B ~ 2,C ~ 4. This underlines 

the importance of the parameter g (and the phase function) for the reflection function calculation. 

Melnikova et al. [7] give the following values for the parameters A,B,C as functions of g: 
A=1.12 + 1.9g, B=198—3g, C=3.514+ 4.8g. (16) 


They were obtained from the parameterizations of results of the solution of the Ambartsumian’s 
nonlinear integral equation at single scattering albedo wp = 1 and the model phase function 


D 

POO F g? F 2g cos 082" 

It follows from Eqs. (16) at g = 0.85: 
Azx4, Ba -—0.57, C759, (18) 


which differ somewhat from our results obtained using p(0), calculated using the electromagnetic 
wave theory [12] for distributions of water droplets. 
We have finally, 
A+ B(E+n)+ Cén + F(0) 
RC, n.p) = i 
A(C +1) 


(17) 


(19) 
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Fig. 2. The relative error of Eq. (19) as a function of the incidence angle at 2 = 0.65 um, aef = 6 um, vı = 15° (a), 
0, =30° (b), 01 = 45° and 60° (c). 


with A,B,C given in Eq. (15) and F(0) = p(0)— p(0). Eq. (19) reduces the calculation of 
R(C,n,p) to the calculation of the phase function p(0) of a given light scattering medium. 


Note that for particular cases (e.g., near-nadir observations) the function F in Eq. (19) can be 
dropped. 


3. The accuracy of approximation 


Let us study the accuracy of Eq. (19) as compared to exact radiative transfer calculations at the 
wavelength 4 = 0.65 um often used in cloud retrieval procedures [13]. The function F(0) can be 


calculated performing the numerical integration in Eq. (5). We, however, applied a different approach 
here. 
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Fig. 3. The same as in Fig. 2 but as a function of the scattering angle 0. 


In particular, we used the fact that the function p(0) can be presented as 


CO 
pO) = X xiP1 (cos 0), (20) 
i=1 
where P;(cos@) is the Legendre polynomial. We see that an infinite discrete set of numbers x; 
completely determines the function p(@). In practice, one can use a finite number M of terms in 
Eq. (20) (M ~ 5Xef, Xef = 274e¢//). Using the addition theorem 
i 
Pi(cos 0) = PROP =n) + Y hiP\(C)PI(H) cosso, (21) 
s=1 


where h; = 2(/ — s)!/(1 +s)! and P(C) are so-called associate Legendre polynomials, we obtain 


BO) = X xP(O)P(—n) (22) 


l=1 
and, therefore, 


[0.0] 
F(0) = X` xı (Pi(cos 8) — (—1)'P(E)P(n)) (23) 
i=1 
where we used the equality: P;(—n) = (—1)'P,(n). 

The function F(0) in Eq. (19) was calculated using Eq. (23) and pre-calculated set of coefficients 
x; obtained with the Mie theory [12] for a given droplet size distribution. 

Relative errors of Eq. (19) with account for Eqs. (15), (23) for spherical water droplets distribution 
at aef = 6 um and A= 0.65 um are presented in Figs. 2a—c. We see in particular that the absolute 
value of the relative error is below 2.5% for most of incidence angles and the observation angle 15°. 
The same is true for smaller observation angles. Therefore, Eq. (19) extends our earlier formulation 
[8] which was applicable only to the case of nadir measurements. 

The errors are larger than 5% only for the case of incidence angles Vo ~ 90° and 3, < 15°, which 
is not of particular interest as far as remote sensing of clouds is concerned. 
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Fig. 4. The same as in Fig. 2 but at aef = 16 um. 


Different points in the vertical direction in Figs. 2a—c represent the azimuth dependence of an 
error. In total, we have calculated 91 azimuth angles [ọ = 0°(2°)180°] for each incidence angle. 
Therefore, Fig. 2 illustrates also the influence of azimuth harmonics on R(¢, n, o ) (for twice and more 
times scattered light). Note that this influence is of increased importance in the rainbow (0 ~ 140°) 
and glory (0 ~ 180°) scattering regions. 

As one might expect, errors increase with the observation angle. However, they are below 5% 
for the observation angles 0-30° (see Fig. 2b) and incidence angles smaller than 60°. They are 
below 7% at incidence angles smaller than 70°. This is also acceptable for most of remote sensing 
applications. Note that errors grow considerably at Jy ~ 90° and V; = 30° (see Fig. 2b). 

For observation angles 45° and 60° only the case of near-nadir illumination (0-30° within accuracy 
5% and 0—40° within accuracy 10%) can be considered in the framework of the approximate formula 
introduced above (see Fig. 2c). 

We present results of calculations given in Fig. 2 as functions of the scattering angle 0 in Fig. 3 
(except for the case 7; =60°). It follows that for large scattering angles (0 > 150°) and observation 
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angles smaller than 45° the accuracy of Eq. (19) is better than 5% (except 0 ~ 180°, where it 
reaches 7%). The same is true till 0 ~ 130° at V, < 30°. 

To check the applicability of Eq. (19) to clouds with larger droplets, we have repeated the same 
calculations at agg = 16 um and give them in Fig. 4. The error is somewhat larger in this case. 
However, it is below 7% at observation angles below 30° and incidence angles smaller than 70°. 
This also indicates the importance of the scattering angle 0, given in Eq. (6) for the accuracy of 
Eq. (19). 


4. Conclusion 


We have proposed approximate Eq. (19) to calculate the reflection function of semi-infinite cloudy 
media. This equation is similar to that presented by Melnikova et al. [7] for the Henyey—Greenstein 
phase function. The accuracy of approximate solution obtained was studied in detail for the water 
droplets distributions f(a) =Aa® exp(—9a/aer) at A= 0.65 um, aef =6 and 16 um. The azimuth was 
varied in the range 0—-180° with step 2° for incidence angles in the range 0—90° and observation 
angles 0-60°. 

It was found that the error of Eq. (19) is below 7% for incidence angles smaller than 70° and 
observation angles smaller than 30°(dep = 6-16 um). For most part of this range the error is below 
5%. Therefore, Eq. (19) can be used in a number of atmospheric optics applications and in particular 
for the satellite cloud retrieval procedures simplification [13,14]. 
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